OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Controlling Malaria Using Livestock-Based Interventions: 
A One Health Approach 

clKkhM updates 

Ana O. M. Gabriela M. Gomes\ Mark Rowland^, Paul G. Coleman^, Clive R. Davies^^ 

1 Instituto Gulbenkian de Ciencia, Oeiras, Portugal, 2 Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, United 
Kingdom 

Abstract 

Where malaria is transmitted by zoophilic vectors, two types of malaria control strategies have been proposed based on 
animals: using livestock to divert vector biting from people (zooprophylaxis) or as baits to attract vectors to insecticide 
sources (insecticide-treated livestock). Opposing findings have been obtained on malaria zooprophylaxis, and despite the 
success of an insecticide-treated livestock trial in Pakistan, where malaria vectors are highly zoophilic, its effectiveness is yet 
to be formally tested in Africa where vectors are more anthropophilic. This study aims to clarify the different effects of 
livestock on malaria and to understand under what circumstances livestock-based interventions could play a role in malaria 
control programmes. This was explored by developing a mathematical model and combining it with data from Pakistan and 
Ethiopia. Consistent with previous work, a zooprophylactic effect of untreated livestock is predicted in two situations: if 
vector population density does not increase with livestock introduction, or if livestock numbers and availability to vectors 
are sufficiently high such that the increase in vector density is counteracted by the diversion of bites from humans to 
animals. Although, as expected, insecticide-treatment of livestock is predicted to be more beneficial in settings with highly 
zoophilic vectors, like South Asia, we find that the intervention could also considerably decrease malaria transmission in 
regions with more anthropophilic vectors, like Anopheles arabiensis in Africa, under specific circumstances: high treatment 
coverage of the livestock population, using a product with stronger or longer lasting insecticidal effect than in the Pakistan 
trial, and with small (ideally null) repellency effect, or if increasing the attractiveness of treated livestock to malaria vectors. 
The results suggest these are the most appropriate conditions for field testing insecticide-treated livestock in an Africa 
region with moderately zoophilic vectors, where this intervention could contribute to the integrated control of malaria and 
livestock diseases. 
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Introduction 

In the last few decades there has been increasing recognition of 
the need for an integrated public health and veterinary approach, 
accounting for the surrounding social-ecological system, to face 
many of the most challenging disease threats: the so-called 'One 
Health' approach [1]. Broadly speaking, animals play an 
important role in the epidemiology of several of the most 
important diseases of man, where they can act as a reservoir 
source for infectious pathogens, and/or a source of blood-meal to 
arthropod vectors of human disease. The recognition of this 
relationship has led to the implementation of human disease 
control strategies targeted at animal populations. These control 
opportunities have been investigated both empirically and 
theoretically. Yet, our knowledge on what determines the public 
health benefits of many of these veterinary interventions remains 
limited. 

A case study of the 'One Health' concept is human malaria in 
regions where its mosquito vectors {Anopheles spp.) also feed on 



animals, since the presence of livestock close to the household can 
affect the rate of vector-human contacts and consequently the risk 
of disease transmission among people. As the Plasmodium malaria 
parasites that infect humans are not infective to livestock, it has 
since long been proposed that animals could be used to divert the 
malaria vector biting from humans, a control intervention known 
as zooprophylaxis [2,3]. However, despite the large number of 
studies performed worldwide for over a century to try to assess the 
value of this strategy in the fight against malaria (reviewed in 
[4,5,6,7,8,9]), the available evidence is still contradictory and no 
consensus exists on the prophylactic effect of animals. Indeed, 
although in several situations the presence of livestock has been 
referred to as a protective factor for malaria vector-human contact 
and/or disease, such as in Papua New Guinea [10,11] and Sri 
Lanka [12], the opposite has been reported in various other 
studies, where livestock were shown to be a risk factor, such as 
Pakistan [13,14], Philippines [15,16], and Ethiopia [17,18] 
(throughout this work the term livestock is used to refer to cattle 
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and other domestic large and small ruminants - buffalos, sheep, 
goats -, as well as donkeys, horses, and swine). 

The apparently contradictory outcomes of the numerous studies 
conducted result from a combination of several possible effects of 
livestock on malaria. On one hand, livestock may divert the blood- 
seeking mosquito vectors from humans, thereby decreasing the 
biting on people [10,11,19] and, as a result, decreasing the 
transmission of the malaria parasite [20] and preventing its 
amplification in people (i.e. the basis for the zooprophylaxis 
concept). But on the other hand, livestock can provide additional 
blood-sources and/or larval breeding sites [21,22,23,24,25], which 
can increase vector survival and/or density [26], consequendy 
increasing the probability of the vector surviving the parasite 
extrinsic incubation period and becoming infectious, as well as 
increasing biting on people [2,6,27]. Additionally, livestock may 
attract more mosquitoes, which, once in the vicinity of the human 
dwellings, may end up biting humans rather than animals 
[14,15,16,18]. The resulting net impact of livestock on malaria 
risk therefore depends on the relative contribution of each of those 
effects. 

In areas where the presence of livestock near people increases 
malaria transmission, an apparently simple solution could be to 
change livestock management in order to deploy the animals away 
from people's houses, between village and vector breeding site [7] . 
However, in Pakistan as well as in some Ethiopian regions, for 
instance, this is not likely to be a feasible strategy, given that 
livestock are such an important source of household income that 
people prefer to keep the animals near their houses to prevent 
them from being stolen [13,28,29,30] and to facilitate husbandry 
practices, such as milking the lactating animals. An alternative 
solution has therefore been proposed: target the non-human host 
of the zoophUic mosquito, by treating livestock with insecticides/ 
acaricides [13] (hereafter referred globally as 'insecticides' for 
simplicity). This strategy has since long been effectively used to 
control ectoparasites and the diseases they transmit to animals 
(and often also to humans), as well as to reduce the direct 
economic losses they cause due to decrease in productivity (e.g. 
lower efficiency of feed conversion, weight gain and milk 
production) [31]. Namely, insecticide treatment of livestock has 
been applied against tsetse flies transmitted animal and human 
trypanosomiasis in sub-Saharan Africa [32,33,34,35], tick-borne 
diseases worldwide (such as anaplasmosis, babesiosis, theileriosis) 
[33,36], and a variety of other biting and/or nuisance flies [37,38], 
mosquitoes [38,39], biting midges [40], mites, and lice. 

The effecti\'eness of insecticide-treated livestock (ITL) against 
malaria was successfully tested by a community-randomised trial 
in Pakistan [41], where the main \'et:tors. An. stephensi and An. 
culicifacies, are highly zoophilic [42]. Notably, following the 
treatment of virtually all domestic animals (93% of the population 
of cattle, sheep and goats) with a solution of the pyrethroid 
deltamethrin applied by sponging, Plasmodium falciparum 
malaria incidence decreased by 56% (95% CI 14%-78%), and 
prevalence decreased by 54% (95% CI 30-69%). Moreover, 
efficacy was comparable to that of traditional indoor insecticide 
spraying but with 80% less costs. Livestock previously infested with 
ectoparasites also improved in weight and milk yield productivity, 
enhancing community uptake of the programme [41]. Additional 
studies have followed to explore whether this strategy could also be 
applied in sub-Saharan Africa, for integrated control of malaria 
and animal trypanosomiasis and tick-borne diseases. Notably, 
bioassays of deltamethrin applied by spot-on and by spray have 
been conducted in Ethiopia [43], and in Tanzania [44], 
respectively, to assess the effects of ITL on the mortality and 
behaviour of malaria vectors. However, despite the encouraging 



results from these bioassays, the impact of ITL on malaria 
transmission at the community level is yet to be formally assessed 
in Africa, where the disease burden is the greatest, but the 
dynamics and determinants of infection differ from Asia. 

A possible concern with ITL is repeUency of mosquitoes, which 
may increase vector feeding on untreated livestock or unprotected 
humans, and make the intervention detrimental. It is known that 
certain insecticides exert not only (1) a toxic or direct insecticidal 
effect, killing mosquitoes that contact with an insecticide-impreg- 
nated surface, but also (2) behavioural avoidance responses. These 
sub-lethal behavioural effects include a) contact-mediated irri- 
tancy, inhibiting mosquitoes from remaining on the treated 
surface, thereby stimulating them to exit prematurely (common 
with pyrethroid insecticides), and b) non-contact or spatial 
repeUency, which acts from a distance of the treated surface 
inhibiting mosquitoes from entering treated areas [45,46]. 
Hereafter, the latter two responses will be referred together as 
repeUency, since any of them could cause mosquitoes diversion to 
another host, in analogy with the shift in host feeding from 
humans to domestic animals that has occasionally been associated 
with the use of pyrethroid-treated nets [47,48,49,50,51]. Addi- 
tionally, a case-control study in the Pokot territory of Kenya and 
Uganda [52] found that people with ITL had a higher risk of 
Visceral Leishmaniasis, suggesting that the insecticide might have 
repeUed sandflies attempting to feed on animals and diverted them 
to feed on humans. Although, to the best of our knowledge, such 
behavioural shift has not been reported for ITL and anopheline 
mosrjuitoes, the possibUity of it occurring should not be 
disregarded and is therefore important to investigate, particularly 
because the most promising insecticides tested on livestock to 
target malaria vectors have been pyrethroids 
[19,39,43,44,53,54,55]. The popularity of pyrethroids is due to 
their high insecticidal action associated with low mammalian 
toxicity [56,57] which makes them safe for both the treated 
animals and for the consumers of animal products. 

An additional concern with using ITL against malaria in Africa 
is that, even in areas where the moderately zoophilic An. 
arabiensis vector (which can easily feed on humans or livestock, 
depending on host abundance and accessibiUty) predominates over 
more anthropophilic vectors such as An. gamhiae s.s., the ITL 
intervention is stUl likely to achieve a smaUer reduction in malaria 
transmission than in Pakistan (and other areas of South Asia), 
where the vectors are highly zoophilic, taking most of their 
bloodmeals upon livestock. A possible way to overcome this 
problem could be to artificially increase the attractiveness of 
insecticide-treated animals to the malaria vector. Although such 
has not been tested in the field yet, the use of synthetic attractants 
to lure anopheUne vectors towards baits or traps and away from 
humans is an area of increasing research [58,59]. 

This work aims to clarify the different effects of livestock on 
malaria and to understand under what circumstances livestock- 
based interventions could play a role in malaria control 
programmes. This was achieved by, firstly, developing a mathe- 
matical model that predicts the apparently contradictory outcomes 
that have been associated with the presence of untreated livestock 
in different ecological settings, and secondly, by expanding the 
model to incorporate insecticide treatment of livestock and fitting 
it to data from Pakistan (where the ITL trial was performed [41]) 
and from Ethiopia (where a field study was conducted [9]) to 
investigate the potential and limitations of ITL. We focus on 
livestock-based interventions, without comparing their effect with 
other malaria control interventions, such as insecticide-treated 
bednets and indoor spraying with residual insecticides. The model 
characterizes situations where livestock by itself can lead to a 
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decrease, increase, or no net impact on malaria transmission to 
humans, and it further indicates that treating livestock with 
insecticide can be a useful complementary tool to control malaria, 
not only in Asia, but also in sub-Saharan Africa. 

Materials and Methods 

Malaria model 

A mathematical model for the transmission dynamics of human 
malaria was developed based on the Ross and Macdonald models 
[60,61], where humans are compartmentalized into either 
susceptible (uninfected and not immune), or infected/infectious 
(SIS model), and mosquito vectors are divided into susceptible 
(uninfected and not immune), exposed/latent (have been infected 
but are not yet infectious) or infectious (SEI model). Here, the 
Ross-Macdonald model is extended by discriminating the feeding 
behaviour of the vector on its alternative hosts: livestock and 
human populations, and by incorporating the treatment of 
livestock with insecticide as a potential novel method to control 
human malaria. The new model explicidy incorporates the effects 
of untreated and insecticide treated livestock on the vector 
population feeding behaviour, mortahty and population density, 
allowing exploration of the impact of livestock-based interventions 
on malaria transmission dynamics. A diagrammatic flow chart of 
the model is presented in Figure 1. Throughout the article, the 
human, vector and livestock populations will be referred to with 
the subscripts h, v and /, respectively. 

The model is formally represented by a .system of ordinary 
difierential equations as follows. For the dynamics of infection in 
the human population, we have 



where = Sf,+If, (total human population). Transmission of 
infection from vectors to humans depends on the number of 
infected vectors per human, IJN^, the vector blood feeding rate 
on any host, a, (the interval between bloodmeals on any host is 1/ 
a), the proportion q of feeds taken on humans (so-called human 
blood index - HBI), the probability h that a human will become 
infected following the bite of an infectious vector, and the number 
of susceptible hosts, S^. Once susceptible humans are infected the 
parasite undergoes a period of latency before infective gametocytes 
appear, but as this period is short compared to the duration of 
infection, it is not represented explicidy in the model [62]. Infected 
individuals, 7^, recover from infection at a rate r, eventually 
becoming fuUy susceptible to re-infection (the average duration of 
infection is 1/r). It is therefore assumed that there is no boosting 
immunity due to repeated infections, as done for simplification in 
earlier zooprophylaxis models [27,63,64,65,66]. Human natural 
mortality and reproductive rates are omitted from the model 
because humans have a long life expectancy relative to other time 
periods used in the model (such as the latent period, infectious 
period and vector life span). We also assume no disease-induced 
death and therefore, the human population size remains constant. 
The disease dynamics in the vector population is represented by 

dS^ ( h , \„ 



where N-^ = Sy+L^+Iy (total vector population). The vector 
dih — ( Q I, fj population comprises only adult female anopheline mosquitoes, 

dt \ N/iJ ' since males do not blood feed. Transmission of infection from 
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Figure 1. Schematic representation of the malaria model, l-iorizontal solid lines denote transitions between epidemiological states, and 
dashed lines represent transmission of infection between human hosts and mosquito vectors. Dotted lines denote vectors feeding on livestock. The 
vector population consists of adult female anopheline mosquitoes. 
doi:l 0.1 371/journal.pone.01 01 699.g001 
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humans to vectors depends on the proportion of infectious 
humans, IiJNh, the vector feeding rate on humans, aq, and the 
probabihty c that a vector will become infected after feeding upon 
an infectious human. Infected latent mosquitoes, Ly, become 
infectious after a sporozoite maturation period (latent period 
= l/(u). Anopheline vectors are assumed to remain infectious 
throughout their life, as usually observed. Infection is assumed to 
have no impact on vector feeding behaviour, reproduction, nor 
mortality, as in most malaria models. Although some effects from 
infection have been described [67], they are not considered in the 
approximation adopted here. 

The vector life expectancy is often about the same order of 
magnitude as the latent period in the vector. Consequentiy, only a 
minority of the infected vector population survives to become 
infectious, and therefore the model must incorporate the class of 
latent vectors as well as vector mortality and recruitment. The 
mortality rate of adult vectors, fi, is assumed to be age 
independent, such that the average vector life-span is We 
consider two implementations for the recruitment rate. One where 
the vector population is kept constant by assuming that 
recruitment and mortality rates are equal p = yU. 

Another where the density-dependent regulation of the adult 
vector population due to competition within the larval stages, 
which depends on the abundance and extent of breeding sites, is 
explicitiy modelled. Following Lord et al. [68] and Kawaguchi et 
al. [66], the recruitment rate of newly emerged female adults 
entering the susceptible class is given by 



/3 = Po-P^^vandp^ 



K' 



K 



where p and are the vector recruitment rate in the absence of 
density-dependence constraints and the strength of the density- 
dependence in recruitment, respectively, and K is the vector- 
carrying capacity of the ecosystem. It is assumed that the number 
and capacity of the breeding sites (and therefore, the vector- 
carrying capacity) remain the same independently of the hosts' 
abundance and availability. For instance, the potential increase in 
breeding sites due to livestock hoof prints is not considered here. 
While density dependence is essential for the systematic investi- 
gation of the zooprophylactic effects of livestock populations with 
dilferent sizes and characteristics, this is no longer a focus on the 
investigation of insecticide treatment later on. As such, variable 
vector population is primarily used in the first part of the Results 
(Untreated Livestock) and the constant vector population is the 
implementation of choice throughout the second part (Insecticide 
Treated Livestock). 

As done for simplification in previous malaria models, we 
assume that vec:tors take one bloodmeal per gonotrophic cycle, 
and therefore, the interval between bloodmeals corresponds to the 
length of the gonotrophic cycle. Similarly, female mosquitoes are 
assumed to feed homogenously with a fixed preference for humans 
and/or animals. 

Formulation of livestock effects 

In the absence of insecticide treatment, the effects of livestock on 
the human blood index, q, follow what has been proposed by Sota 
and Mogi [63] and are defined by: 



NhAh 



which can be simplified to: 



1 + 



Ni_Ar 

NhAh 



where Ah and A] are the proportional availabilities of the human 
and livestock hosts, respectively, and can take any value between 0 
and 1, inclusive. The term availability encompasses all the factors 
that can influence the likelihood of the vector feeding on a given 
type of host, when two types of alternative hosts are present in 
equal numbers. Namely, these factors include the accessibility of 
each host to the vector (which can vary with distance between 
vector breeding sites and location of humans/livestock at night, 
whether located indoors or outdoors, under a bednet or not, or 
hvestock enclosed inside a shed or not), and on the intrinsic 
propensity to feed upon humans versus animals (anthropophily 
versus zoophily), and to feed in the location where the host resides 
(endophagy versus exophagy), which can be modified by vector 
genetics and learning. In the presence of insecticide treatment, the 
expression for the human blood index is generalized as 



1 

NiAt{\- 



(3) 



Nh Ah 



where e is the proportion of livestock population treated with 
insecticide, hereafter referred as treatment coverage, and a is the 
diversion probability, defin(;d as the probability that a host-seeking 
mosquito will be diverted away from (a > 0, repellency) or towards 
(a<0, attractancy) an insecticide-treated animal. Therefore, the 
insecticide treatment of livestock only affects the human blood 
index if the intervention has some diversion effect upon the 
vectors, either repellency or attractancy (see Text Sl.l for more 
details). 

The baseline mortality rate of the vector is decomposed as being 
the sum of the minimum mortality rate (^i^) due to causes other 
than searching for a bloodmeal host (i.e. mortality due to hazards 
during the act of feeding on a host, the gestation period, the search 
for oviposition sites, and the underlying aging process), and the 
mortality due to searching for a bloodmeal host {p.^. The search- 
related mortality is assumed to be proportional to the length of the 
searching period, which is inversely related to the abundance and 
availability of potential blood meal hosts. These assumptions 
follow previous models by Saul [27] and Killeen and Smith [65]. 

When no livestock are treated with insecticide, the expression 
for the vector mortaUty rate therefore becomes 



1 



{NhAh+N,A,)j 



a. 



where the last term is the search-related mortality, p^. Parameter^' 
is a factor to scale the proportional availabilities [Ah, Ai) of hosts to 
the mosquito vectors into absolute availability values. As in most 
previous malaria models, it is assumed that the feeding success of 
malaria vectors is independent of the density of vectors per 
available host [69]. 



PLOS ONE I www.plosone.org 



4 



July 2014 I Volume 9 | Issue 7 | e101699 



Controlling Malaria with Livestock Interventions 



When livestock are treated with insecticide the mortality rate is 
generalized as 

J« = /<»,+J«S + J«k 

, / 1 \ f e(l-oc)N,A! \ (4) 

In the term, as seen for the human blood index, if there is 
repellency (a>0) it is as if the availability of livestock became 
reduced by the proportion ea, which corresponds to the 
proportion of bites attempted on a given animal that will be 
diverted to another animal or human host. This wiU cause an 
increase in the time it takes for the vector to find a bloodmeal host, 
with consequent increase in the search-related vector mortality. 
Conversely, if there is attractancy (a < 0), it is as if the availability 
of livestock became increased by the proportion ea, which 
corresponds to the proportion of bites attempted on a given 
animal that were diverted from another animal or human host. 
This will decrease the time for the vector to find a bloodmeal host, 
thereby decreasing the search-related vector mortality. The fi/^ 
term accounts for the direct lethal effect of insecticide applied on 
livestock, and is a function of the vector biting rate on livestock, the 
treatment coverage, e, the diversion probability, a, and the 
insecticidal probability, k. The daily biting rate, a, needs to be 
included in the expressions for jx^ and /i^-, since the additional 
mortalities, either due to searching for a bloodmeal host or due to 
attempting to feed on insecticide-treated livestock, are only 
suffered by the vector when it attempts to blood feed (see Text 
SI. 2 for more details). Our model assumes that the insecticide 
effects (diversion and insecticidal probabilities) are constant, 
therefore reflecting average values of what would be observed 
throughout the year. 

Simulations 

The system of equations (l)-(2) was analysed symbolically for 
the derivation of endemic equilibrium solutions (see Text S2) and 
numerically for the simulation of dynamical trajectories over time. 
Numerical integration was performed using BERKELEY MA- 
DONNA V. 8.3.9, with the built-in method /oitrt/i order Runge- 
KuUa. The equilibrium solutions were further explored with 
MATLAB V. R2011a. 

We first investigate the effects of untreated livestock in malaria 
transmission and then move to explore the impact of treating 
livestock with an insecticide that has lethal and possible 
diversionary effects (repellency or attractancy) upon malaria 
vectors. For this purpose, a range of simulations was performed 
with system (l)-(2), focusing on scenarios of endemic Plasmodium 
falciparum malaria. 

Threshold derivation 

We also determined the threshold conditions required for 
persistence of malaria, by analyzing the equilibria of the model 
represented by system (l)-(2). The average number of secondary 
cases generated by a single infectious individual introduced in a 
population of fully susceptible individuals, is known as the basic 
reproduction number, denoted by Rq [61,70]. This threshold 
quantity expresses the transmission potential of an infectious 
disease and must exceed unity for the infection to be maintained in 
the population. The expression for Rq was derived by linearization 
around the disease-free equilibrium (DEE), based on the next- 
generation operator approach [71,72]. We then explored the 
impact of ITL on Rq for different intervention scenarios. Namely, 



by setting Rq-\, we obtained the critical proportion of the 
Uvestock population that must be treated with insecticide, and 
assessed how this critical coverage would be affected by the 
insecticide diversionary properties. 

Parameterization 

Parameters values for the untreated livestock model were 
obtained directly or derived from the literature and are provided 

in Table 1. The effects of insecticide treatment were explored 
using parameter values that were either extracted or derived from 
empirical data from the index studies in tiu; North- West Frontier 
Province of Pakistan (ITL trial conducted by Rowland et al. [41]) 
and in the Konso district of South- West Ethiopia (field study by 
Franco [9]), or from previous studies within or near the area of the 
index studies, as listed in Table 2. See Text S3 for details on 
parameterization. 

Results 

Untreated livestock 

Here we explore the effects that varying the abundance and/ or 
availability untreated livestock could have on different outcome 
measures of malaria transmission. All the simulations used 
parameter values as listed in Table 1, unless otherwise specified. 

Figure 2 shows the vector population density and prevalence of 
human infection over time as li\ estock are introduced in a setting 
where previously only humans and no livestock were present. 
Simulations were performed assuming a fix human density 
(A';, = 100) and 1 head of livestock per person [Bffi — 0.25). The 
proportional availability of livestock to vectors was the same as 
that of humans for all plots in this figure (4; = 0.5), illustrating the 
case of a moderately zoophilic vector, like An. arabiensis in 
Ethiopia. Additional scenarios of host density and availability were 
also explored. 

Firsdy, we simulate a modified model with the best case scenario 
where the vector population is kept constant (A'„(t) = 1000) by 
assuming that recruitment and mortality rates are the same {p = fl) 
(black line), and secondly, the carrying capacity was set to a higher 
level (^ = 5,000 to 100,000) and the vector population density 
increased from its initial equilibrium (A'„(0) = 1 000) towards 
carrying capacity (coloured lines). As we would expect, in the 
case of constant vector density, the introduction of livestock leads 
to consistent reductions in the prevalence of human cases by 
diverting vector feeds to livestock (black line). Overall, the higher 
the numbers and/ or availability of the introduced livestock, the 
stronger is the predicted zooprophylactic effect on malaria 
transmission. When the vector population density is allowed to 
increase, however, the prevalence of human cases might increase 
(coloured lines). For a given density and availability of livestock, 
the higher the carrv'ing capacit}' is in relation to the initial \'cct()r 
population density, the higher the vector density and consc'quently 
malaria transmission levels in the new endemic equilibrium. In all 
simulated scenarios in Figure 2 the system reaches a new 
equilibrium in less than 3 years after the introduction of livestock. 

Figure 3 examines various outcome measures of malaria 
transmission that characterize the new endemic equilibrium that is 
reached under a range of relative livestock to human density {9ni 
varying from 0 to 1), when the proportional availability of livestock 
to vectors is either the same (^4^ = 0.5) or nine times higher 
{Al = 0.9) as that of humans, the latter resembling a scenario of a 
highly zoophilic vector, ~ An. culicifacies in Pakistan. The 
outcome measures investigated include the human blood index 
(HBI, designated as q in our model), daily overall vector mortality 
(fi), vector density (A^„),daily entomological inoculation rate (EIR), 
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Table 1. Paranneter values for modelling the effects of untreated livestock on malaria. 



Symbol 


Definition 


Value 


[Reference] 


a 


Vector daily biting rate on any host 


0.5 


[89] 


b 


Probability that humans become infected from the bite of an infectious vector 


0.04 


[90] 


c 


Probability that vectors become infected after biting on an infectious human 


0.3 


[90] 


r 


Human daily recovery rate from infection (1/average duration of infection) 


0.05 


[29,30,90] 


m 


Daily rate at which infected mosquitoes become infectious (1 /latent period} 


0.07 


[91] 


/I 


Overall average vector dally mortality rate (/(,„ + /fv) 


Varied 


Derived 




Vector daily mortality rate In absence of available livestock 


0.1 


[89] 


Mm 


Vector dally minimum mortality rate when there are no hazards due to search for a bloodmeal host 


0.05** 


[89] 


J"» 


Vector daily mortality rate due to searching for a bloodmeal host 


Varied** 


Derived 


p 


Overall average vector dally recruitment rate 


Varied 


Derived* 


Po 


Vector daily recruitment rate in the absence of density-dependence constraints 


Varied 


Derived* 


Ps 


Strength of the density-dependence in recruitment (/day) 


Varied 


Derived* 


K 


Carrying capacity of the vector population (/ha) 


10^ to 10^ 


_ 


NM 


Initial vector density, prior to change In livestock abundance and/or availability {/ha) 


10^ 




Wft 


Human density (/ha) 


100 




em 


Relative density of livestock:humans (A///A//,) 


0 to 20 




A, 


Proportional availability of livestock to vectors 


0 to 1 




Ah 


Proportional availability of humans to vectors { = ^-A^ 


0 to 1 




q 


Proportion of vector bloodmeals on humans (Human Blood Index) 


0 to 1 


Derived 


j 


Scaling factor to transform proportional availabilities Into absolute availabilities 


Varied 


Derived 



*For simulations with constant vector population density: p = ^; for variable vector density: p = Po(l —N^/K) and = f.if,K /{K — Nv(0}). 
**The relative magnitudes of fi^ and were varied in a sensitivity analysis. 
doi:l 0.1 371 /journal.pone.Ol 01 699.t001 



and prevalence of infection in liunians (//,). Tlie EIR is the number 
of infective mosquito bites received by a human per unit time, 
estimated multiplying the daily human-biting rate (HBR) by the 
proportion of mosquitoes with sporozoites in their salivary glands 
(I^/N^). The HBR is the total number of mosquito bites received 
by a human, per day, and is calculated as the product of the 
number of vectors per human and the number of daily bites on 
humans per vector (HBR = (A'j,/A'/,)flHBI).The figure illustrates 
the effects of livestock on decreasing the human blood index while 
decreasing vector mortality (Figure 3A,B) and increasing vector 
population density (Figure 3C,D). The combination of these 
effects may lead to situations where the presence of livestock 
increases, decreases, or has no significant impact on malaria 
transmission (all other panels in Figure 3). The introduction of 
livestock is predicted to have a zooprophylactic effect, i.e. decrease 
malaria transmission, in two situations. One is that the vector 
population density does not increase as a result of livestock 
introduction. The other is that although the vector population 
density increases as a result of livestock introduction, the livestock 
numbers and availability to vectors are sufficiently high, such that 
the increase in vector density is counteracted by the diversion of 
bites from humans to animals (Figure 3). Otherwise, the 
introduction of livestock is predicted to increase malaria transmis- 
sion. 

Impact of vector search-related mortality on the effects of 
untreated livestock. For the purpose of illustrating the model 
behaviour, the simulations for untreated livestock assume that the 
vector-search mortality when no livestock are available (i.e. when 
N[ = 0 or Ai = 0), has the same value as the vector minimum 
mortality rate (/i, = /(^ = /i/,0.5 = 0.05/day). A sensitivity analysis 
was done to explore the impact of different relative magnitudes of 



the vector search-related mortality (Figure SI in Text S4.1). If 

the vector search-related mortality is already negligible before 
livestock are introduced, then introducing livestock will have no 
impact on the vector mortality, and will simply decrease HBI, 
consequently decreasing malaria transmission. Conversely, if the 
vector search-related mortality is considerable, introducing live- 
stock can considerably decrease vector mortality, which can 
increase the proportion of vectors surviving the extrinsic incuba- 
tion period to become infectious, and thereby counteracting the 
decrease of the HBI due to diversion of mosquito bites from 
humans to livestock, consequently increasing malaria transmission. 
After a certain threshold of livestock density, further increasing 
their abundance produces negligible reduction on vector mortal- 
ity. 

Insecticide treated livestock 

To explore the effects of ITL on malaria the model was fitted to 
P. falciparum malaria transmitted by the highly zoophilic An. 
culicifacies in Pakistan and the more anthropophUic An. arabiensis 
in Ethiopia. Parameter values are listed in Table 2. The main 
differences in the malaria transmission parameters between the 
Asian and African settings are as follows. In Ethiopia, livestock 
were 8. 1 times more abundant, although with an estimated 56.8 
times lower availability to the main malaria vector, than in 
Pakistan, resulting in a predicted HBI over 4 times higher in the 
African than in the Asian setting. Additionally, the estimated 
duration of the latent period in vectors was slightly shorter, while 
the vector life expectancy was 75% higher in Ethiopia than in 
Pakistan. The initial density of vectors per human and the 
probability of infection in vectors were set to be, respectively, 3.3 
and 13.6 times higher in Pakistan than in Ethiopia. 
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Table 2. Paranneter values for modelling the effects of insecticide-treated livestock on malaria. 





Symbol 


Definition 


Value 




[Reference] 








Pakistan 


Ethiopia 


Pakistan 


Ethiopia 


a 


Vector daily biting rate on any host (1/gonotrophic cycle) 


0.4 


0.4 


[92] 


[93,94] 


b 


Probability that humans become infected from the bite of 
an infectious vector 


0.5 


OS 


[95,96] 


[95,96] 


c 


Probability that vectors become infected after biting on 
an infectious human 


0.95 


0.07 






r 


Human daily recovery rate from infection 
(1 /average duration of infection) 


0.05 


0.05 


m* 


*** 


0) 


Daily rate at which infected mosquitoes become infectious 
(1 /latent period) 


0.057 


0.064 


Derived from 
[9,91] 


Derived from 
[9,91] 


H 


Overall average vector dally mortality rate + + 


Varied 


Varied 


Derived 


Derived 




Vector daily natural mortality rate in the absence of ITL 
(1 /natural life expectancy) 


0.22 


0.12 


Derived from 
[41,92] 


Derived from 
[93,94,97] 




Vector daily minimum mortality rate when there are no 
hazards due to search for a bloodmeal host 
(1/vector maximum life expectancy) 


0.11**** 


0.06**** 


_ 


_ 




Vector daily mortality due to searching for a bloodmeal host* 


0.11**** 


0.06**** 


Derived 


Derived 




Vector daily mortality due to the direct lethal effect of 
insecticide applied on livestock 


Varied 


Varied 


Derived 


Derived 


P 


Overall average vector daily recruitment rate 


= 


= f^ 






Ny 


Vector density (/ha) 


5000 


1500 






Nh 


Human density {/ha) 


100 


100 


- 


- 


Bm 


Relative density of livestock:humans (W/A//,) 


0.14 


1.13 


[41] 


[9] 


Oai 


Relative availability of livestock:humans (/\///\/,) 


53.24 


0.938 


Derived from 
[98] 


Derived from 
[99] 


A, 


Proportional availability of livestock to vectors (Oai /{^ + Oai)) 


0.982 


0.484 


Derived 


Derived 


Ah 


Proportional availability of humans to vectors 


0.018 


0.516 


Derived 


Derived 


1 


Proportion of vector bloodmeals on humans* 


0.118 


0.485 


Derived 


Derived 


J 


Scaling factor to transform proportional availabilities into 
absolute availabilities 


Varied 


Varied 


Derived 


Derived 


e 


Treatment coverage: proportion of livestock population 
that is treated with insecticide 


Oto 1 


0 to 1 






k 


Insecticidal probability 


0.1 (0 to 0.9) 


0.1 (0 to 0.9) 


Derived from 
[54] 




a 


Diversion probability (a>0, repellency; a<0, attractancy) 


0 to 1 


-1 to 1 






Malaria vectors: An. culicifacies in Pakistan, An. arabiensis in Ethiopia. 

*Parameter values pre-intervention that will be affected if livestock are treated with an insecticide with diversion properties. 
**Values chosen to produce malaria prevalence similar to the observed In the index study areas. 



***M. Rowland unpublished data. 
****The relative magnitudes of /i^, and //„, were varied in a sensitivity analysis. 
doi:l 0.1 371 /journal.pone.Ol 01 699.t002 



It was assumed that, prior to the ITL intervention, an endemic 
equilibrium of malaria transmission had been reached and, as in 
most previous zooprophylaxis models [27,65,66], vector popula- 
tion density was at its equilibrium level, and remained constant 
throughout the intervention (i.e. vector recruitment and mortality 
rates are the same). We therefore consider the scenario where the 
insecticide has no impact on the overall vector population density. 
Thus, the beneficial impact of an ITL intervention with a non- 
diversionary insecticide is assumed to be due only to the decrease 
on vector survival caused by the toxic insecticidal effect, and 
consequent reduction in the proportion of vectors that become 
infectious. When the insecticide additionally has some repellent 
properties, there is some beneficial effect from increasing the 
vector search-related mortality, which partially counteracts the 
increase in vector bloodmeals on humans. Conversely, when there 



is attractancy, there is the greater benefit of decreasing the 
bloodmeals in humans, which counteracts the decrease in vector 
search mortality. 

Impact on malaria prevalence. We started by exploring 
the predicted impact of ITL on the prevalence of human infection. 
This is represented in terms of the prevalence ratio (PR), which is 
defined as the ratio between the prevalence under a given 
coverage of insecticide-treated livestock (e) and the prevalence pre- 
intervention. The proportional reduction on the pre-intervention 
prevalence is given by 1 -(prevalence ratio). 

Simulations were initially performed to estimate the coverage of 
treated livestock (e) and insecticidal probability {k) required to 
obtain the 54% reduction in P. falciparum prevalence observed in 
the Pakistan ITL trial, for the Pakistan and the Ethiopian 
simulated scenarios, assuming the use of an insecticide with no 
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diversion properties (a = 0, Figure 4). For any given intervention 
effort, ITL is predicted to cause a stronger reduction in malaria 
prevalence in Pakistan than in Ethiopia. Nevertheless, the same 
reduction in prevalence could be achieved in Ethiopia, if using 
higher treatment coverage and/ or a product with stronger or 
longer lasting insecticidal properties. For instance, for the scenario 
oik = 0.1 (estimated value from Pakistan data, as detailed in Text 
S3. 2), the predicted treatment coverage required to obtain the 
observed reduction in prevalence (PR = 0.46) is e=15% in 
Pakistan and 25% in Ethiopia (Figure 4). 

We also investigated whether by increasing the attractiveness of 
insecticide treated livestock to vectors it would be possible to 
obtain in Ethiopia the same reduction in prevalence as observed in 
the Pakistan trial, with the same interx^ention effort used in the 
Asian setting (Figure 5). To achieve in Ethiopia the same 
PR = 0.46 with similar coverage as predicted ft)r Pakistan 
(£= 15%, for A = 0.1, assuming no repellency) would require an 
attractancy of 20%, while with attractancy of 10% or 30%, the 
coverage would be approximately 19% or 13%, respectively. 

We then explored how the coverage would be affected if the 
insecticide had a repellency effect upon vectors, and what might be 
the repellency probability above which the intervention could 
become deleterious, by causing prevalence to increase above the 
pre-intervention level. Not surprisingly, the intervention benefits 
considerably decrease if the insecticide has repellency properties. 
Considering again the case of the estimated A; = 0. 1 (Figure 5), to 
achieve the observed reduction in prevalence (PR = 0.46) with the 
93% coverage that was actually applied in the Pakistan trial, the 
model suggests that a repellency probability of ~ 1 7 % would need 
to be acting in Pakistan. If that same repellency level was acting in 
Ethiopia, the required coverage was predicted to be 60%. For 
repellency above 1 7 % in Pakistan or above 2 1 % in Ethiopia, the 
achieved reduction in prevalence is expected to be always smaller 
than the observed (i.e. the prevalence ratio, PR, would always be 
>0.46), even if all livestock are treated (e= 1). The intervention 
would become deleterious (PR >1) for repellency above 20% in 
Pakistan and above 28% in Ethiopia (Figure 5). The smaller the 
c<)\'cragc (for a gi\Tn k), or the greater the /{ (for a given coverage), 
the higher is the repeUence threshold above which ITL will start 
becoming detrimental (PR >1) (Figure S2 in Text S4.2). 

Threshold phenomena. The derived basic reproduction 
number for the malaria model is given by: 



Nr (aq) he CO 
Ni, rn (£u + /i) 



(5) 



where q is given by expression (3) and ji is given by expression (4). 

By setting /?q = 1 in (5) we see that the critical proportion of the 
Uvestock population that must be treated with insecticide in order 
to interrupt malaria transmission is 



To explore how repellency (a > 0) or attractancy (a < 0) would 
impact Ec numerical simulations were performed (Figure 6). The 

stronger the insecticidal probability (^), the smaller is the critical 
proportion of treated livestock (e^) required to potentially reduce 
Rq below unity, for any given a in Ethiopia, and for a<0.4 in 
Pakistan. In the Asian scenario, for a>0.4, above a certain 
treatment coverage and insecticidal probability, there could be a 
shift from Rq< 1 to Ro> 1 . For instance, for a = 0.4 and k = 0.4, Rq 
becomes less than 1 if coverage is above 36% and below 90%, 
while for coverage above 90% then Rg increases to greater than 1, 
Similarly, for a = 0.5 and k = 0.6, Rq is reduced to less than 1 for 
coverage between 29% and 70%, but for coverage above 70% the 
Ro becomes above 1 . For any given k, the stronger the repeUence, 
the higher is the critical coverage, while the stronger the 
attractancy, the lower is the critical coverage. Furthermore, for 
any given k, with or without repellency, the critical coverage is 
always higher for Ethiopia than for Pakistan (Figure 6). 

Impact of vector search-related mortality on the effects of 
insecticide-treated livestock. The baseline simulations for 
ITL assume that the background vector search-related mortality 
(pre-livestock treatment) has the same value as the vector 
minimum mortality rate (/i, = /J„, = /Jq.0.5). Additional simulations 
were done to explore the sensitivity of the findings to alternative 
search-related vector mortahty values (Figure 5 and Figure 6 
can be contrasted with Figure S3 and Figure S4 in Text 
S4.2.1, respectively). 

Although there is uncertainty about its exact value, the relative 
magnitude of the background vector search-related mortality will 
only affect the intervention impact if the insecticide has 
diversionary properties. Namely, decreases in the background 
search mortahty will counteract the only benefit of repeUence 
(which was an increase on the search-associated vector mortality), 
and consequently decrease the beneficial impact of an ITL 
intervention. In general, the smaller the background search-related 
mortality, the stronger is the detrimental effect of any given 
repellency probabUity (a>0) on malaria prevalence or Rq, and 
conserjuently, the greater is the coverage rerjuired to achieve a 
given reduction in prevalence or Rq, and the lower is the 
repeUence threshold above which the intervention would become 
deleterious (and vice-versa). For instance, comparing the baseline 
scenario with the worst-case scenario of nuU background vector 
search-related mortality, the repeUence threshold would decrease 
from 20% to 13% in Pakistan and from 28% to 19% in Ethiopia 
(Figure S3 in Text S4). This relationship becomes however 
increasingly non-linear with increase in the insecticidal effect (k) of 
a treatment with repellency, namely in Pakistan (Figure S4 in 
Text S4). For a given attractancy probabUity (a < 0), the smaller 
the background vector-search related mortality, the stronger are 
the intervention benefits, and consequently, the smaller is the 
coverage required to achieve a given reduction in prevalence or Rq 
(Figure S3 and Figure S4 in Text S4). 



no\ 4^(aq'fbc+rm ] -r[ tu + 2( /(,„ + 



2ra{\-q')k 



.(6) 



where is the HBI in the absence of insecticide treatment: 

This expression for Sc is valid for the best-case scenario 

regarding repeUence and the worst case scenario regarding 
attractancy, i.e. when using an insecticide without any diversion- 
ary effects upon vectors (a = 0). 



Discussion 

By combining a mathematical model with field data we have 
explored the different effects that livestock can have on human 
malaria in areas where the disease is transmitted by zoophUic 
vectors, allowing us to understand under which circumstances 
fivestock-based interventions could play a role in malaria control 
programmes. 

Our model predicts that the presence of untreated livestock wUl 
have a zooprophylactic effect in two scenarios. One is when vector 
population density does not increase as a result of livestock 
introduction. The other is when although the vector population 
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Figure 2. Temporal effect of introducing livestock in a setting 
with endemic malaria. Effect of introducing livestock in a setting 
where only humans were present, when: A/„ remains constant (black 
line), and when increases until reaching a maximum, which depends 
on the carrying capacity, /C (increasing from /C= 5,000 (green) to 100,000 
(red)). Nv,(0)= 1000, Nh= 100 and /\/ = 0.5: the availability of livestock to 
vectors is the same as that of humans; Oa>/ = 0.25 (1 head of livestock 
per 4 persons). To achieve the same initial equilibrium A/„ (and /f,) 
for various K values, the vector recruitment rate in the absence of 
density-dependence constraints was set to vary accordingly: 
PD = ^i,K/{K — Nt,(0)). Other parameters are as in Table 1. 
doi:l 0.1 371/journal.pone.Ol 01 699.g002 

density increases, the numbers and availability of livestock to 
vectors are sufficiently high (such that the resulting diversion of 
bites from humans to livestock can counteract the increase in 
vector density), or the vector mortality related with host-search 
pre-livestock introduction was sufficient low (such that introducing 
livestock causes no significant decrease on the already small search 
mortality). Otherwise, the introduction of livestock is predicted to 
increase malaria transmission. 

These results are in agreement with the insights from two 
previous zooprophylaxis models [27,63]. Namely, Sota & Mogi 
[63] also identified as key determinants of the beneficial versus 
detrimental effect of untreated livestock on malaria transmission 
whether the vector population had reached its maximum possible 
density prior to livestock introduction, and whether the density 
and/or availability of animal hosts were sufificientiy high. It is 
worthwhile mentioning that these features are captured by both 
the present model and the Sota & Mogi [63] model although the 



two works differ in the approach used to model the potential 
detrimental impact of livestock on malaria transmission. The 
present work explicitiy models the effect of animal or human hosts' 
abundance and availability on vector mortality, with consequent 
impact on the dynamics and density of adult vectors. Instead, Sota 
& Mogi [63] assumed a constant vector mortality rate, and 
modelled the effect of hosts abundance and availability on the 
probability of successful blood feeding of the vector, with 
consequent impact on the number of eggs laid and density of 
adult vectors in the future generations. Aside from the work by 
Sota & Mogi [63], two other previous zooprophylaxis models that 
addressed the effect of untreated livestock on malaria transmission 
[27,65], have also explicitiy modelled the effect of animal or 
human hosts abundance and availability on vector mortality. Saul 
[27] also highlighted that the effect of untreated livestock on 
malaria greatly depended on the magnitude of the search-related 
vector mortality: when this is significant, increase in livestock 
density could lead to increased malaria transmission, which is 
consistent with our results. 

Regarding the insecticide-treatment of livestock, when using an 
insecticide without diversionary properties, any given intervention 
effort is predicted to achieve a stronger reduction in malaria 
transmission in a setting with highly zoophHic vectors (exemplified 
by Pakistan) than with the more anthropophilic An. arabiensis 
(illustrated by Ethiopia), as expected. Yet, the same reduction in 
malaria prevalence could be achieved in Ethiopia, if treating a 
high proportion of the livestock population with a product that has 
stronger and/ or lost lasting insecticidal effect than what was used 
in Pakistan. The predicted intervention effort required to achieve a 
given reduction in prevalence with a non-repellent insecticide, is 
however, surprisingly low, and most likely unrealistic. In the 
Pakistan trial, a 54% reduction in prevalence was obtained, 
following treatment of 93% of the livestock population in the trial 
villages (cattle, goats, and sheep) [41]. Our results suggest that, to 
achieve the observed reduction in prevalence with such high 
treatment coverage, the insecticidal effect would need to be 
extremely small. 

When accounting for a possible repellency effect of the 
insecticide, the expected benefits of the intervention decrease 
considerably in both settings, requiring more realistic parameter 
values to obtain the results observed in the Pakistan trial. 
Repellency threshold probabilities were identified above which 
the intervention could become detrimental, increasing the 
prevalence of human infection above the pre-intervention levels. 
For repellency probability below those thresholds any vector 
diversion to humans was predicted to be overcompensated by the 
insecticidal (direct lethal) effect and the increased search-related 
mortality of the mosquitos attempting to blood feed on insecticide- 
treated animals. Within that range of repellency probability for 
which ITL is likely to stiU reduce malaria prevalence, a greater 
benefit may be observed in Pakistan or in Ethiopia, depending on 
the repellency and coverage levels. 

The results indicate that repellency has a stronger detrimental 
impact on malaria (prevalence or in Pakistan than in Ethiopia, 
and therefore, it would take a smaller level of repellency for ITL to 
start becoming deleterious in the Asian setting. Above the 
repellency threshold the intervention becomes always more 
detrimental in settings with higher availability of livestock to 
vectors, like in Pakistan and other settings with highly zoophilic 
vectors. 

The repellency level of the insecticide applied to animals can 
thus have an important effect on the intervention outcome. For a 
given treatment coverage, the stronger and/or longer lasting the 
insecticidal effect, the higher is the repellency threshold above 
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Figure 3. Effect of altering the relative livestock to human density, for different vector density scenarios, at the new endemic 
equilibrium. Comparing a scenario where the availability of livestock to vectors is the same as that of humans (left, /A/ = 0.5) versus where it is 9 times 
higher than that of humans (right, /A/ = 0.9). Along the x-axis, representing OMi = N/Nh, the livestock density N/ is varied relative to a fixed human 
density N/, = 100. N„(0)=1000. Effect of introducing livestock when: remains constant (black line), and when Wi, increases until reaching a 
maximum, which depends on the carrying capacity, K (coloured lines: K increasing from 5,000 (green line) to 100,000 (red line)). The effects of 
introducing livestock on the human blood index (HBI) and on the vector mortality rate (/i) are independent from the vector density scenarios (A, B). 
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The vertical line in the left panels highlights the new endemic equilibrium that is reached after the introduction of 1 head of livestock per 4 persons 
(Ojv/ = 0.25), corresponding to the end of the timeline in Figure 2. Other parameters are as in Table 1. 
doi:1 0.1 371 /journal.pone.01 01 699.g003 

transmission, obtaining field estimates for this component of vector 
mortality is an important challenge that future research should 
address. 

The repellency threshold above which the intervention might 
become detrimental could be as low as 13% in Pakistan and 19% 
in Ethiopia, if assuming aU livestock population is treated with an 
average direct insecticidal effect of 10%, under the worst case 
scenario of null vector search-related mortality. The smaller the 
treatment coverage, and/or the stronger the insecticidal effect or 
the search-related mortality, then the higher the repellency level at 
which ITL can still be safely used. 

To our knowledge, this is the first modelling approach that 
explicidy explores the potential effects of repellency and attrac- 
tancy in the context of ITL and malaria transmission. The present 
work is an improvement in relation to previous malaria models of 
the impact of applying insecticide on animals [27], on animal 
sheds [66], or on bednets [65]. None of the former two models 
[27,66], explored a repellent or attractant effect of the insecticide, 
and although work by KUleen and Smith [65] has looked at 
repellency and livestock applied to an African setting, it did so in 
the context of insecticide-treated bednets and diversion of malaria 
vectors to humans and/or untreated cattie, without referring to 
insecticide-treatment of catde. 

Considerations on modelling repellency 

The present work assumes that when a mosquito tries to bite on 
an insecticide-treated animal and is repelled, it will be diverted to 
bite on another host. Nonetheless, it could be that the mosquito is 
not able to find a successful bloodmeal and does not feed in that 
night, ending up either feeding only on the foUowing night, or 
dying earUer. The impact of repellency on vector mortality is 

Pakistan Ethiopia 




Coverage (e) Coverage (e) 

Figure 4. Predicted impact of Insecticide Treatment of Livestock on malaria prevalence, without diversion {a = 0). This figure shows the 
combination of values of coverage and insecticidal probability required to achieve a given prevalence ratio (PR: prevalence with ITL / baseline 
prevalence). Blue line: PR = 0.46 (like the observed in the Pakistan trial); White line: PR = 0; Dashed line: k = OA, as estimated for the Pakistan trial. The 
colour bar shows the scale of PR values, from 0 to 1. Other parameters are as in Table 2. 
doi:1 0.1 371/journal.pone.01 01 699.g004 



which the intervention starts becoming detrimental. Additionally, 
when considering doing ITL interventions with high treatment 
coverage of the livestock population, researchers should be aware 
that the higher the intervention coverage, the greater the 
detrimental effect from a given repellency level, and the greater 
the benefits from reducing repellency. A small decrease in 
repellency could gready improve the intervention benefits, with 
the effect being greater in scenarios with more zoophUic 
mosquitoes. Interestingly, this is the opposite from the case of 
insecticide-treated nets (ITNs), where repellency could be bene- 
ficial in some circumstances. Namely, the greater the proportion of 
the human population covered with ITNs the greater are the 
expected benefits from repellency, and, conversely, the smaller the 
coverage the greater the likelihood that malaria vectors might be 
diverted from ITN-protected people to those unprotected (if the 
density and/ or availability of animal hosts to the mosquito vector 
are small) [65]. 

In general therefore, the smaller the repellency, the greater the 
benefits of an ITL intervention. The benefits in settings with 
moderately zoophUic vectors, such as An. a7'ahiensis in sub- 
Saharan Africa, could be further improved by artificially 
increasing the attractiveness of livestock to the malaria vector. 

If the insecticide has diversionary properties upon the malaria 
vectors, the magnitude of the vector mortality related with host 
searching was predicted to considerable affect the model results. 
Namely, the smaller the vector search-related mortality pre- 
intervention, the stronger are the insecticide diversionary effects 
upon malaria prevalence or R^, be it the detrimental effect of a 
given repellency probability on transmission, or the reduction in 
transmission obtained with a given attractancy probability. Given 
the influential role of the vector search-related mortality upon the 
effects of untreated and insecticide-treated livestock on malaria 
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Figure 5. Predicted impact of Insecticide Treatment of Livestock on malaria prevalence - with repellency {y.>0) or attractancy (s<<0) 
for Ar= 0.1. This figure shows how the diversionary properties of the insecticide affect the coverage required to achieve a given prevalence ratio (PR: 
prevalence with ITL / baseline prevalence). Blue line: PR = 0.46 (like the observed in the Pakistan trial); White line: PR = 0; Red line: PR = 1 (above which 
treating livestock increases malaria prevalence). Along the y axis, a is varying from no diversion (o; = 0) to maximum repellency (a=l) or maximum 
attractancy (a= — 1). The colour bar shows the scale of PR values, from 0 to =1 1 in Pakistan and up to ==5 in Ethiopia. Other parameters are as in 
Table 2. 
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captured by the model, since repellency reduces the availability of 
treated livestock, increasing the time required to find a bloodmeal 
host and consequently increasing the vector search-related 
mortality. The impact of repellency increasing the interval 
between bloodmeals is something that could be explored by 
extending the model to explicitly account for that possibility. 

We also assume that the probability of vectors being repelled to 
humans after attempting to bite on livestock, depends only on the 
proportion of livestock population that is treated with insecticide 
(coverage,£), on the repellency probability of the insecticide (a>0), 
and on the relative number and availability of livestock or human 
hosts. Additionally, the model assumes that repellency and 
coverage are independent. In reality, however, the occurrence of 
a repellent effect can depend on additional factors such as 
characteristics of the: a) insecticide (chemical compound, formu- 
lation and concentration); b) intervention (concentration of the 
insecticide on the animal's coat, which will eventually decrease 
with time after application); and c) mosquito vector [73]. Also, the 
insecticide concentration is likely to be heterogeneous throughout 
the animal's surface, and the place where the mosquitoes land on 
the animals can therefore be determinant. 



With regards to the mode of action of insecticides applied on 
livestock depending on the properties of the insecticide itself, some 
pyrethroids are more toxic to vectors than repellent (e.g. 
deltamethrin, used in the Pakistan ITL trial), other pyrethroids 
are more repellent than toxic (e.g. permethrin), and other classes of 
insecticides (e.g. organophosphates) are just toxic and non- 
repeUent. Yet, even the typical toxic deltamethrin tends to be 
repellent at low dosages. Namely, as the applied dose of 
deltamethrin decays over time it goes from being toxic to non- 
toxic but repellent and then to just repellent. 

Due to this, a big concern during the Pakistan ITL trial [41] was 
that mosquitoes would be repelled onto humans as the dosage of 
deltamethrin decayed, but it appears malaria was still controlled 
because the insecticide was reapplied regularly before there was 
too much decay. This explanation is consistent with the findings 
from the present work where, on one hand, when accounting for 
repellency the model results are more compatible with the 
observed Pakistan trial results, than when assuming that the 
insecticide had no diversion effect. On the other hand, the 
predictions suggest that the stronger and/or longer lasting the 
insecticidal effect, the highest is the repellency threshold above 
which the intervention is likely to become detrimental. Addition- 
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Critical coverage {ej Critical coverage (e^) 

Figure 6. Critical proportion of ITL as a function of the insecticidal [kt, and diversionnary effect (a). The lines show the combination of 
values of coverage and insecticidal probability required to achieve Ro= 1, above which Rq will be decreased below 1, for a given diversion probability 
(a). Black line: s( = 0, no repellency or attractancy (is the same as the white line in Figure 4); Red lines: a>0, repellency increasing from 0.1 to 0.5 (top), 
at intervals of 0.1; Green lines: a<0, attractancy increasing from —0.1 to -0.5 (bottom), at intervals of 0.1. Other parameters are as in baseline 
simulations (Table 2). 
doi:1 0.1 371 /journal.pone.01 01 699.g006 



ally, at the high treatment coverage applied in the trial, there is 
only a small difference between the repellency level with which the 
observed reduction in prevalence would be achieved, and the 
repellency threshold. This supports the hypothesis that the 
intervention effort applied in the Pakistan trial was sufficiently 
high to make the repellency effects non evident. 

Increasing livestock attractiveness to vectors 

By increasing the attractiveness of insecticide treated animals to 
malaria vectors, it could be possible to ffirther enhance the impact 
of ITL in malaria control in settings with more opportunistic 
vectors, such as An. arabiensis in Ethiopia, as shown in this work. 
This could eventually enable extending the geographic regions 
where ITL might reduce malaria burden, to include also areas 
with more anthropophilic vectors, such as An. gambiae s.s., the 
most competent malaria vector in sub-Saharan Africa. Given the 
potential benefits this could bring, it would be worthwhile further 
exploring this hypothesis in future work. 

In practice, however, insecticides tend to be non-attractant (i.e. 
neutral or repellent). Therefore, to artificially increase livestock 
attractiveness would require developing an insecticide that has also 
attractancy properties (in addition to its toxic insecticidal effect), or 
alternatively, treat livestock with an attractant substance on top of 
applying a standard insecticide. Although this may sound 
somewhat speculative, it is not much different from what has 
been successfully tested in other systems, where synthetic 
attractants have been applied to baits or traps to increase their 
attractiveness to tsetse flies [74,75,76], anopheline mosquitoes 
[58,59], and other insects of medical and veterinary importance 
[77]. 

Regarding possible detrimental implications of artificially 
attracting more mosquitoes into livestock, these are likely to be 
minimal. Attracting a mosquito to a cow does not necessarily 
mean the mosquito will succeed in biting/blood feeding as it may 
be killed or knocked down by exposure to the insecticide before 
taking up blood, and that is usually the case, namely with 



pyrethroids. The expected reduction in mosquito survival due to 
increased exposure to the insecticide toxicity should actually lead 
to less biting. Therefore, it is unlikely there would be additional 
disease burden or economic costs, as long as the attractancy would 
be specific for malaria vectors and would not cause increased 
number of biting ffies or other arthropods that are vectors of 
pathogens to livestock, and would also not cause a reduction in the 
animal's blood through excessive biting that could decrease milk 
or meat yield. 

Optimizing insecticide-treated livestock interventions 

It is important to highlight that, although we explored the 
impact that treating livestock with insecticides could have on 
malaria transmission, this intervention has been traditionally used 
with a veterinary purpose, to control tsetse flies, ticks and other 
ectoparasites, and the diseases they transmit to animals, improving 
livestock health and productivity, such as milk and meat yield. 
Therefore, when evaluating the cost-effectiveness of the interven- 
tion, both the animal health benefits and the public health benefits 
need to be captured ('One health'). Given the potential double side 
benefits of veterinary interventions like this, and given the central 
role of livestock in poor tropical settings, to control human disease 
and improve livestock health will have disproportionate economic 
impact that needs to be captured, as accounting for it could 
promote the wider implementation of the intervention. Namely, if 
the costs of ITL are allocated to the human health and the animal 
health sectors in proportion to the benefits, the intervention might 
be profitable and cost-effective for both sectors. Here lies a 
challenge to the Public Health community, which will require 
strengthened collaboration with the Animal Health community. 

In addition to the animal health/ productivity benefits, ITL uses 
much less insecticide than traditional malaria control methods, 
such as indoor spraying of houses with residual insecticide, making 
ITL very cheap from a human disease control perspective. In the 
Pakistan trial, sponging livestock with deltamethrin was shown to 
achieve a reduction in malaria burden similar to indoor residual 
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spraying but with l^>Q% less campaign costs. Furthermore, if 
accounting for the increase in milk production by the treated 
cattle, associated with clearance of tick infestations, the economic 
gain would he enough to cover all insecticide and labor costs [41]. 

The mathematical model developed here could be used to 
examine the economic aspects of the 'One Health' approach to 
disease control, encompassing both himian and animal health 
benefits at a societal level. The model provides a framework for 
quantifying the benefits of ITL as a reduction in the human health 
burden (expr(-sscd as prevented DALYs, Disabilit)' Adjusted Life- 
Years), associated reduction in health care costs (expressed as $), as 
weU as the improvements in animal health and productivity 
(expressed as $). The cost-effectiveness of ITL, accounting for both 
the human and veterinary benefits, could then be compared with 
other interventions that deliver only human health benefits, such 
as indoor residual spraying (IRS) and insecticide-treated bednets 
(ITNs), and the relative attractiveness of ITL across epidemiolog- 
ical settings and animal production systems examined. 

The use of any animal-based intervention for malaria control 
win only be a component of the broader integrated malaria control 
approach, and will have to be deployed alongside case detection, 
treatment and prevention. The relative importance of animal- 
based interventions within the broader approach will vary between 
settings. Also, the adoption of any recommended intervention is 
intimately related to the socio-economics of the setting and it is 
therefore vital to understand the drivers for adoption by the target 
population. 

In Pakistan very high treatment coverage was achieved with a 
free campaign and the animal owners were enthusiastic because 
they could see the benefits of tick elimination and improved mUk 
and meat yield [41]. Previously to the campaign the insecticide 
treatment of livestock for ectoparasites was normally ad hoc done 
by householder according to perceived need, which would lead to 
only partial coverage at any one time. Therefore, a subsidized 
campaign approach is recommended, similarly to the externally 
funded campaigns of IRS. 

Although the empirical evidence for Africa is lacking, some 
inferences can be made from tsetse control work. In particular, the 
experience from controlling human sleeping sickness in southeast 
Uganda by targeting the cattle reservoir of the human infective 
parasite shows that large scale campaigns can also reach very high 
(>80%) coverage levels with insecticide and trypanocidal treat- 
ment. Additionally, reducing the volume of insecticide, and so the 
price of ITL treatment, through restricted apphcation protocols to 
target insecticide use to those areas of the cattle where tsetse or 
anopheline mosquitoes preferentially feed, have the potential to 
drive routine ITL adoption by small-holder farmers [35,43]. To 
drive the private uptake of ITL usage, farmers need to see a direct 
benefit to their animals. Experience from the sleeping sickness 
work shows that for effective control through ITL it is important 
that the insecticide products used work against both ticks and the 
human disease vectors (such as synthetic pyrethroids), as tick 
control is often the main motivation for farmers to use ITL [34]. 

ITL is particularly useful for malaria control where vectors, in 
addition to bloodfeeding on livestock, are (or have became) 
exophagic (feeding outdoors, therefore escaping to ITNs exposure) 
and/ or exophilic (resting outdoors, and thereby evading IRS). 

One cannot rule out that long term and intensive use of ITL 
may lead to selection for anthropophily, with a consequent shift in 
preference from animals to humans (assuming that host preference 
is determined by genetic polymorphisms [78,79,80]). Therefore, 
changes in the HBI (as a proxy for host preference) should be 
monitored in regions where repeated campaigns are undertaken 
[14,41]. Additionally, selection for anthropophily could be 



countered by combining ITL with indoor strategies to control 
anthropophiUc and endophUic mosquitoes, like ITNs and IRS 
[19]. 

At the time of the field studies in the settings to which to which 
the ITL model was parameterized (Konso region of Ethiopia and 
NWFP in Pakistan) most people were not using bednets. Future 
work could expand the present model to investigate the use of 
livestock-based interventions alongside ITNs or IRS, to provide 
additional insights to the potential impact that combining these 
strategies might have on malaria transmission. 

A concern inherent to any vector control intervention based on 
insecticides is the potential development of resistance. Namely, 
pyrethroid resistance is becoming increasingly wide spread across 
anopheline mosquitoes [81,82] and several other arthropods that 
feed on livestock, such as ticks [83,84]. It has been argued that the 
treatment of livestock with pyrethroids is not likely to induce 
stronger selection pressure for resistance in malaria vectors than 
insecticide-treated bednets or indoor residual spraying of houses 
and cattie sheds, but nevertheless, appropriate monitoring of the 
vector populations is required if wide scale and long term ITL 
interventions are implemented [14,41,54]. 

It has also been recommended that research efforts should 
target the identification of alternative non-pyrethroid insecticides 
for livestock treatment [54]. Possible candidates have recentiy 
been suggested from the avermectins class of insecticides, which 
have since long been used in veterinary and human medicine 
against several helminths and arthropod pests [85] and were latest 
shown to be also toxic to anopheline mosquitoes. Namely, feeding 
on bovine blood treated with ivermectin reduced survivorship and 
fecundity oi An. gambiae s.s. and^n. arahiensis [86,87] and may 
possibly also inhibit the sporogony of P. falciparum as it was 
recentiy shown in treated humans [88]. Another promising 
avermectin is the more recent eprinomectin which has similar 
antihelminthic and ectoparasiticidal action as ivermectin in cattle, 
but with much less mammary excretion, allowing its use in 
pregnant and lactating animals, on the contrary of ivermectin 
[87]. Any of these avermectines could overcome the problems of 
pyrethroid resistance as well as repellency upon malaria vectors, 
and could be administered as part of mass livestock vaccination 
campaigns, simultaneously benefiting animal and human popula- 
tions. Additionally, while pyrethroids can only be administered 
topically, both ivermectin and eprinomectine are available 
topically (as pour-on) and also as injectable formulation (subcu- 
taneous administration), which could surmount the difiiculty faced 
with pyrethroids of achieving high enough concentrations of 
product throughout the animal's skin. Malaria vectors would 
however need to bite the animal and take a bloodmeal to be 
exposed to the insecticide, but every biting mosquito would be 
exposed and die more promptiy, therefore requiring a smaller 
dose, compared to pyrethroids. Further studies are needed to 
assess the effects of livestock treated with the recommended dose of 
ivermectin or eprinomectine upon wild populations of malaria 
vectors. 

Conclusions 

A mathematical model was developed to predict the different 
effects of untreated and insecticide-treated livestock in malaria 
outcomes in different regions. Similarly to previous work, our 
model indicates that the zooprophylactic effect of untreated 
hvestock depends on whether 1) the pre-existing malaria vector 
population had reached its maximum density, 2) livestock 
abundance and availability to the vector is suflficientiy high, and 
3) vector mortality related with host-search pre-livestock intro- 
duction was sufEcientiy low. We additional find that, as expected, 
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the insecticide-treatment of livestock is likely to be more beneficial 
to humans in settings with highly zoophilic malaria vectors as in 
Pakistan and other areas of South Asia, than in settings with 
moderately zoophilic vectors, as An. arahiensis in sub-Saharan 
African. Nevertheless, the intervention could also substantially 
decrease malaria burden in the latter settings, under certain 
conditions, as illustrated here with the predictions for Ethiopia. 
Namely, in regions with moderately zoophilic vectors the benefits 
of the intervention will be maximized if 1) treating most of the 
livestock population with a product that has a stronger or longer 
lasting toxic insecticidal effect than what was used in the Pakistan 
trial, and that has htde (ideally nuU) repellency effect (such as the 
non-pyrethroids ivermectin or eprinomectin), or 2) if the 
attractiveness of the treated animals to malaria vectors could be 
increased. 

It is hoped that this work may lead to increasing awareness 
about the non-Knear effects of livestock on malaria transmission, 
and to the implementation of a community-based trial of 
insecticide-treated livestock in an African region where An. 
arahiensis predominates, and where this strategy could potentially 
contribute to the integrated control of human malaria and 
livestock diseases. 
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